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Abstract 

To clarify the relationship between phase-stable acceleration (PSA) and unlimited relativistic 
acceleration (URA) (Phys. Rev. Lett. 104, 135003 (2010)), an analytical relativistic model is 
proposed in the interactions of the ultra-intense laser and nanometer foils, based on hydrodynamic 
equations. The dependence of the ion momentum on time is consistent with the previous results 
and checked by PIC simulations. Depending on the initial ion momentum, relativistic RPA contains 
two acceleration processes: phase-stable limited relativistic acceleration (PS-LRA) and URA. In 
PS-LRA, the potential is a deep well trapping the ions. The ion front, i.e., the bottom, separates 
it into two parts: the left half region is PSA region; the right half region is PSD region, where the 
ions climb up and are decelerated to return back. In PS-LRA, the maximum ion energy is limited. 
If the initial ion momentum large enough, ions will experience a potential downhill and drop into 
a bottomless abyss, which is called phase-lock-like position. URA is not phase-stable any more. 
At the phase-lock-like position the ions can obtain unlimited energy gain and the ion density is 
non-zero. You cannot have both PSA and URA. 
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Laser-ion acceleration has been an international research focus l|-|4|, however it is still 
a challenge to obtain mono-energetic proton beams larger than lOOMeV. As a promis- 
ing method to generate relativistic mono-energetic protons, radiation pressure acceleration 
(RPA) has attracted more attention 3 61 and becomes dominant in the interaction of the 

n 

ultra- intense laser pulse with nanometer foils. The phase-stable acceleration [7|] predicted 
that the energy spread can be improved in the interactions of nanometer-foils with circular- 
polarized laser pulses. With thin-shell modeljs], Bulanov and coworkersjs] pointed out that 
the ions can obtain unlimited energy gain by RPA in the relativistic limit. Yan and cowork- 
ers tried to predict the ion energy distribution with a self-similar hydrodynamic theory {9)]. 
However, it is nonrelativistic and under the plasma approximation which allows V«i?^0 
and rii — n e = satisfy together, where E is the acceleration field and rii (or n e ) is the ion (or 
electron) density. However, can the ions obtain unlimited energy gain in the phase-stable re- 
gion or is the phase-stable acceleration still possible in the relativistic limit? How about the 
relationship between them or what are the critical conditions for them? To give clear answers 
of the questions, an analytic self-consistent relativistic fluid model is proposed to describe 
the relativistic radiation pressure acceleration and to recheck the unlimited ion-acceleration 
region. 

The ion acceleration in the interaction of ultra-intense laser and nanometer foils contains 
two stages: the hole-boring process and the radiation pressure acceleration. Here the tran- 
sition of them is assumed steady and the instability of the acceleration sheath is suppressed 
well, which can be realized for specially designed targets 10l. 111]. In the hole-boring process, 
the ion velocity can reach M/i&Q, [l3], which is the hole-boring velocity and also the initial 
velocity of the ions in the radiation pressure acceleration. The initial time, to is when the 
compression layer is detached from the foil. It is decided by the target thickness and Uhb- 
With the initial conditions: Uhb and to, the dependence of the ion energy on time can be 
obtained and consists with the results of thin-shell model|3l and has been checked by PIC 

n 

simulations [3]. Depending on the laser intensity, the initial ion momentum will determine 
two different acceleration processes: the phase-stable limited relativistic acceleration (PS- 
LRA) and the unlimited relativistic acceleration (URA). When the initial ion momentum 
is smaller than the critical one, the potential is a deep well trapping the ions. The accel- 
eration mode is PS-LRA and contains the phase-stable acceleration region (PSA) and the 
phase-stable deceleration region (PSD). The well is separated into two half- regions: the left 



half- region and the right half- region by the bottom, which is the ion front. The left half- 
region is PSA region, where the electric field is positive and ions coast down to the bottom. 
While the right half-region is PSD region, where the electric field is negative and ions climb 
up the potential uphill and are decelerated to return to the bottom. The deceleration is 
also phase-stable. No matter in PSA or PSD region, the maximum ion energy is limited 
and ascertained by an analytical formulation. Since PSA and PSD are separated by the 
ion front where the ion density is zero, the ions in two regions can not exchange from each 
other. If the initial ion momentum is large enough, the ions can get across the potential 
uphill and experience a potential downhill and then drop into the bottomless abyss which 
is the phase-lock-like position. If the ions can reach the phase-lock-like position, they can 
obtain unlimited energy gain as Bulanov and coworkers pointed out [8]. The acceleration 
mode is URA and not phase-stable any more. The electron density is smaller than the ion 
density and the electron front increases with time. The phase-lock-like position of URA is 
the limiting ion front and the ion density is non-zero as time tends to zero. You cannot 
have both PSA and URA in the same acceleration process. Therefore, the maximum ion 
energy is finite in the phase-stable region and URA is not phase-stable any more. However, 
no matter in PS-LRA or URA, the plasma tends to neutral as time tends to infinite. 

For convenience, the physical parameters: the time, t, the ion position, x, the ion velocity, 
v, the electron field, E, the electric potential, <p, the plasma density, n, and the light speed, 
c, are normalized as follows: r = ut,x = xk,u = v/c,E = E/Eo,<f> = (p/{po,h = n/n,Q, 
where n represents rij (or n e ), is the reference density, u is the light frequency, k = u/c 
is the wave number, c is the light speed, E = k<p , e(p = •y em m e c 2 and 7 em is the maximum 
electron energy. Here e is the elemental charge. 

With reference to the results given by Mako and Tajima in Ref. I], in the self-similar 
state, the density distribution of ions is assumed as: 

n k = ^(l + i>) a ,k = l,...,N, (1) 

where the subscribe & stands for the ion species, Qk is the charge number of the kth species 
ion, the index a depends on the laser intensity and the target thickness and discussed in the 
Ref. Q. 

With the transformation: £ = x/r, the normalized continuity and motion equation of 



ions and Poisson's equation are given as: 



d\nn k du k 
[u k - i) ^ = — — , 



-p (EQ k h k - n e 



where /3 k = < ^ kl ^ me ; M k is the mass of certain ions, 
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Solving Eq. (jSJ), the ion velocity satisfies: 
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Mfc,0 = U k + ^ — (x - Xo) 
la 



(4) 



and the potential in the ion region is given by: 



(x~X0- 2au kfll k fi) 

1= 4aJ k 



where 7 fc , = (1 - u 2 k0 ) 1/2 , xo = x(«fc,o)» 

X = ll du k , 
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and u k: o is the hole-boring velocity given by [12|, Il3| : 
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where a 2 = 0.732/ 10 i8 W / cm 2A 2 tm , Iioisw/cm 2 is the laser intensity in unit of 10 18 W/cm 2 , and 
= (1 — (uhb/c) 2 )^ 1 / 2 . The beginning time To is when the compressed ion and electron 



layer is detached from the foil and given by 
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at £ = 0. 

In order to obtain the dependence of the ion velocity on time, it is need to solve Eq. 
From Eq. fll]), the following differential equations of two variables are given: 

( d Pk 2aV k {l+p 2 fl 2 
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dt 



t \3ap k V k ^l+p 2 - (2a + l)t 

dVk = -2aV k 

dt 3ap k V k ^l + p 2 - {2a + l)t 
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FIG. 1. (Color online) Comparison of our analytical fluid model with thin-shell model and Esirke- 
pov et al.'s simulations for a/a ~ 0.1, a = 316, d = A = lfj,m, uq = 49n c = 5.5 x 10 22 /cm 3 , and 

a = l^r 

where pk = UkJk is the normalized ion momentum and Vk satisfies: 

Vk = f P* dt _ t Vk . (10) 

Using matlab function odell3 to solve Eq. (JH]) with the initial time Tn and initial velocity, 



Mfc,o, the dependence of the ion energy on time has been calculated [14(. As an example, 



Figure [T] shows the comparison of our analytical fluid model with the thin-shell model and 
PIC simulations [3|. The results of our model consist well with that of the PIC simulations. 
When t is larger than 40 times of the laser cycle, our results are a little larger than that 
of the PIC simulations. One of the reason is the loss of laser energy and the decreasing of 
the electron temperature for large r in the simulations, while the electron temperature is 
assumed to be a constant in our model. 
Combing Eqs. (TJQ) and (j5J), it is obtained: 

1 (x~Xo- 2au k,oltfl^ 

With Eqs. ([3]), @ and (JS]), the electron density in the ion region is written as: 

n e = (l + 0ir + -^^, (12) 
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Equation (11 2p shows the plasma can not be quasi-neutral at a finite time. However, the 
plasma tends to neutral as the time tends to infinite. It is obtained: 



lim n e = SQ fc ra fc , (13) 
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With Eq. (fTTj) and rifc(£ = £jj) = 0, the possible maximum ion energy at the ion front 
satisfies: 



Pk.m 3 fPkfi 
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Ik 2d Pk= / 7 fc 2 ^Pfc + 2aPfc,o7fc, > ( 14 ) 



where 7| = 1 + pi, ptm = «fc,m/ ~ u tm * s ^ ne normalized limiting momentum of ions at 
the ion front: f = £ i)f and p M = u kfi / Jl - u 2 kfi . 

Different from the nonrelativistic case, it contains two acceleration modes and depends 
on the initial conditions: Ukfl and a in the relativistic case. The critical condition is decided 
by: 

rPk,o __3 1 r+00 

/ 7 fc 5 dPfc + 2ap M 7fc = / i 1 + ^Y m dx « 2.622. (15) 
Jo Jo 
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First, J Pfc '° 7 fc 2 <ipfc + 2apfc j0 7fc < 2.622 for $ 0.5064 and a = 2, which requires 
a < 203 for no = 49n c , d = A = Ifim. Figure |2] shows the dependence of the plasma density, 
electric field and potential on £ for n = 10 22 /cm 3 , a = 2 and a = 70. In this case, it 
is divided into two regions: the phase-stable acceleration region for < £ < y and the 
phase-stable deceleration region < £ < 1. 

(I) In the phase-stable acceleration region, < £ < the electric field > and 
the electron density is larger than the ion density. The maximum ion momentum pk, m i is 
limited and given Eq. ([14"]) at the ion front ascertained by Eq. with = Uk >rn i, 
where Uk, m i = Pk,mi/ + v\ m \- At ^ ne i° n fr° n t, the ion density is zero. The difference of 
the electron density and ion density decreases with time and tends to zero. 

The potential shown by Figure [2] (b) in < £ < £j j gives an intuitionistic explanation 
of PSA. The ions coast down the slope of the potential, and the gradient, i.e., the electric 
field, becomes gently as the ions come to the bottom of the potential although few can reach 
there. Therefore the ions at higher potential will obtain more acceleration and the energy 
spread is improved. That is PSA. Different from the real gliding process, the ions can not 
pass through the bottom and climb up since the ion front is the limiting point and the ion 
density tends to zero at the bottom. 



(II) In the phase-stable deceleration region, £ it f < £ < 1, the eletric field E < and the 
electron density is larger than the ion density too as shown in Figure [2j In this region, the 
ion momentum p k satisfies: p k}Tnl < p k < p k , m2 , where p kjm2 = p k (£, = 1) and given by Eq. 
01]) with £ = 1. All the ions in this region are decelerated to £ = The absolute value of 
the electric field decreases with the decreasing £ and is zero at the ion front. Therefore the 
deceleration is also phase-stable. 

With Figure 121(b), for ions, the potential in < £ < 1 is a mountain with height of about 
35GeV which is far larger than the maximum ion energy of about \Jpl n ^ + 1 — 1 ~ 7GeV 
at £ = 1. Therefore, the ions will be decelerated at the potential uphill. Since the gradient, 
i.e., the value of electric field increases with the potential height, the deceleration is also 
phase-stable. As point out above, the limiting point is still the ion front £ = and the 
ions in PSD region can also not pass through the bottom into PSA region. The ions in PSA 
and PSD region can not exchange from each other because of the zero-density dividing point 

In this case, the maximum ion momentum is p k)in i at the ion front £ij < 1. Therefore, it 
is called phase-stable limited relativistic acceleration (PS-LRA). The ions in the two phase- 
stable regions can not exchange from each other and the ion front £jj is the dividing line. 
Combing the condition of PS-LRA and the above discussion about PSA and PSD, the ions 
with an initial momentum of p k ^ not large enough to get across the potential at £ = 1 will 
drop in PSA region or PSD region and obtain a finite maximum energy ascertained by Eq. 

If the initial ion momentum satisfies: 

fPkfi _3 1 

/ 7fc 'dp k + 2op M 7j? > 2 - 622 - ( 16 ) 
Jo 

the ions will not experience a potential well in PS-LRA and will coast down from the potential 
slope and drop into the bottomless abyss at £ = 1 as shown by Figure |3] (d). £ = 1 is called 
phase-lock-like position and the ions can obtain unlimited energy gain as shown in Figure 
|3](b). It is called unlimited relativistic acceleration (URA), which is not phase-stable any 
more. Eq. (1TB]) satisfies for p kt0 > 0.5064, i.e., a > 203 from Eq. ([7]) for a — 2, no = 49n c = 
5.5 x 10 22 /cm 3 and d = A = lum. Figure |3] shows the dependence of the plasma density, 
electric field and potential on £ for a = 316, d = X = lum, no = 49n c = 5.5 x 10 22 /cm 3 , and 
a = 2. In this case, the electron density is smaller than the ion density and the acceleration 
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FIG. 2. (Color online) Phase-stable limited relativistic acceleration (PS-LRA) contains two regions: 
phase-stable acceleration region for < £ < £jj ~ 0.87 and phase-stable deceleration region for 
0.87 < £ < 1. (a)The density of ions and electrons for different time and the ion momentum VS 
the self-similar variable £. At the ion front £jj = 0.87, the ion density is zero. (b)The potential 
VS f. (c) and (d) The electric field E VS £. E > for < £ < and E < for &j < £ < 1. 
Here, n = 10 22 /cm 3 , a = 2 and a = 70. 

is not phase-stable any more. In all the region, the electric field increases with £ and is larger 
than zero. At a finite time, the electron front, where n e = 0, £ e j < 1, which is 0.971, 0.994, 
0.99997 at r = 5, 20, 2000 separately. Therefore the possible maximum the ion momentum 
is shown by Figure [3] (a) and (b). 

As discussed above, due to the initial ion momentum pkfl large enough, the ions can reach 
URA region. Therefore the ion momentum and the field have a sharp increase and tend to 
infinite, the ion density becomes a non-zero constant in £ £ (1 — 5, 1] as r — > oo, where 5 is 
an infinitesimal. As shown in Figure [3] (b), at the phase- lock-like position and the limiting 
ion front £ = 1, the ion can obtain unlimited energy gain and the ion density is non-zero. 
It is similar with the unlimited phase-lock ion acceleration as pointed out by Bulanov and 
coworkers [8] in the relativistic limit. In URA region, it is found that (I) the unlimited ion 
acceleration requires the initial ion momentum is large enough and should meet Equation 
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FIG. 3. (Color online) Unlimited relativistic acceleration (URA) with the phase-lock-like position 
£ = l.(a)The density of ions and electrons for different time and the ion momentum VS £. The 
ion density is no-zero at the limiting ion front: £ = 1. The electron density is smaller than that of 
ions at any finite time. (b)The enlargement of (a) near £ = 1. £ = 1 is the phase-lock-like position 
where the ion momentum tends to infinity, (c) and (d) The electric field E and the potential <j) 
VS £. E > for all the region and the acceleration is not phase-stable any more. Here a = 316, 
d = A = 1/J,m, no = 49n c = 5.5 x 10 22 /cm 3 , and a = 2. 

( I16p : (II) it is not phase-stable any more; (III) the phase-lock-like position is £ = 1, the 
limiting ion front; (IV) the ion density at the limiting ion front is non-zero. 

In the conclusion, it has been given an analytical relativistic fluid model to describe the 
relativistic radiation pressure acceleration with the initial parameters from the hole-boring 
stage. The dependence of the ion velocity on the acceleration time can be obtained and is 
consistent with that of thin-shell model and PIC simulations. There are two acceleration 
modes: PS-LRA and URA with a critical initial ion momentum ascertained by an explicit 
formulation. In PS-LRA, the ions are trapped in a deep potential well and the maximum 
ion energy is limited and the ion front is the well bottom and < 1. URA is not phase- 
stable any more and there is a phase- lock-like position in it. At the phase- lock-like position, 
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corresponding to the relativistic limit, the ions can obtain unlimited energy gain and the 
ion density is non-zero as time tends infinite. Although the unlimited ion acceleration can 
not be reached at any finite time, the ions can be accelerated to any large energy if the 
laser pulse is long enough. As an important result, you cannot obtain both PSA and URA. 
Therefore, if the laser parameters are large enough to obtain URA, the energy spread must 
be lost. If one wants to improve the energy spread with PSA, the maximum ion energy is 
limited. 
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